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We study the quantum localization phenomena for a random matrix model belonging 
to the Gaussian orthogonal ensemble (GOE). An oscillating external field is applied on 
the system. After the transient time evolution, energy is saturated to various values 
depending on the frequencies. We investigate the frequency dependence of the satu- 
\ rated energy. This dependence cannot be explained by a naive picture of successive 

independent Landau-Zener transitions at avoided level crossing points. The effect of 
quantum interference is essential. We define the number of Floquet states which have 
large overlap with the initial state, and calculate its frequency dependence. The num- 
U: ber of Floquet states shows approximately linear dependence on the frequency, when 

i the frequency is small. Comparing the localization length in Floquet states and that 

in energy states from the viewpoint of the Anderson localization, we conclude that the 
Landau-Zener picture works for the local transition processes between levels. 
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The diffusive phenomena have been actively studied in quantum systems with a complex level 
structure in the presence of a periodic external field. Under a periodic perturbation, when a system 
shows so-called quantum chaotic nature, it is widely believed that energy diffusion of the system 
eventually ceases and the energy saturates at a finite value. So far, only limited number of these 
systems have been actually studied. The kicked rotator model^-* and the kicked top modeP^ are 
typical quantum models which show the quantum chaotic nature under a periodic external field. 
The difference between classical and quantum energy diffusion has been discussed by investigating 
these models. In the quantum kicked rotator model, starting initially from the ground state, the 
system evolves diffusively only for a finite time. After then the diffusive time-evolution ceases and 
the average energy is saturated to a finite value. That is, the system cannot absorb energy after 
that time. This is a remarkable difference from the classical diffusion, where the energy increases 
unlimitedly. The mechanism of this quantum localization in the model was discussed in the context 
of the Anderson localization.^'^) Recently the kicked rotator model is experimentally realized using 
Hydrogen and Sodium atoms, and the quantum localization is observed. ^"^^ 

It is known that random matrices well describe characteristics of the level statistics of the 
* machida@spin.t.u-tokyo. ac.jp 
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systems which show chaotic nature in the corresponding classical models. Therefore we study 
the quantum localization process of random matrices to extract universal features. In this paper, 
we investigate the time evolution of a system whose Hamiltonian is taken from the Gaussian or- 
thogonal ensemble (GOE) with a periodic perturbation, and try to capture the universal features 
on the relation between the saturated energy and the frequency of the perturbation. So far, several 
types of time-dependent random matrices are adopted to study the quantum diffusion process. ^^"^^^ 
Wilkinson calculated the energy diffusion constant using the Landau-Zener transition formula^^"^^) 
when the perturbation acts on the system linearly in time.^^^ He clarified the difference of diffusion 
constants between GOE and the Gaussian unitary ensemble (GUE). Wilkinson neglected the effect 
of interference of quantum phase, which plays a crucial role in the quantum localization in the 
presence of periodic perturbation.^^) Cohen and Kottos considered the energy diffusion before the 
saturation in the Wigner's banded random matrix (WBRM) model with an oscillating perturba- 
tion, and calculated the diffusion coefficient as a function of the amplitude and frequency of the 
perturbation.^^) Taking these previous studies into account, we study how the Landau-Zener pic- 
ture is related to the quantum localization by investigating the dependence of the saturated energy 
on the frequency of the perturbation. 

We first numerically demonstrate that the quantum localization occurs in the system, which 
cannot be explained by a naive application of an independent Landau-Zener picture to the system. 
In order to understand the quantum localization intuitively, we consider the overlap between Flo- 
quet eigenstates and the initial ground state, and find that the overlap exponentially decays in the 
Hilbert space spanned by Flouqet eigenstates, which directly means that the ground state is local- 
ized in this representation. We define the relevant number of Floquet eigenstates which has a large 
overlap with the initial ground state. This number was originally introduced as a quantity which 
corresponds to the Lyapunov exponent in the quantum kicked top model by Haake et al?^ The 
dependence of the relevant number on the frequency is numerically investigated. We numerically 
find that this number linearly depends on the frequency of the perturbation in small frequency 
regime. Employing a plausible proportional relation between the localization length in the energy 
space and that in the Floquet space, which is also adopted in the case of the kicked rotator model, 
and taking account of this linear dependence, we find the Landau-Zener mechanism works in the 
view of the Anderson localization. 

This paper is organized in the following way. In §2, we explain the random matrix model and 
the numerical method of time evolution. In §3, we study the frequency dependence of the saturated 
average energy. In §4, we discuss the quantum localization by introducing the number of relevant 
Floquet states. Finally in §5, we give summary and discussion. 
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We consider time-reversal invariant systems under an oscillating external field. The random 
matrix ensemble appropriate for describing the spectral statistics of these systems is the Gaussian 
orthogonal ensemble (GOE). Random matrices well describe characteristics of energy spectra of 
complex qTiantum systems which have no conserved quantities. 

We shall consider the total Hamiltonian given by 

H{t) = Ho + X{t)V, (1) 

where Ho denotes the non-perturbed part of the Hamiltonian, and X{t)V is the perturbation part 
with a time-dependent parameter X{t). We take Ho and V from the ensemble of the GOE random 
matrices with the dimension N. This model corresponds to the realistic systems such as billiard 
systems with a time-dependent boundary, or complex spin systems with a time-dependent external 
field. Matrix elements Hoij and Vij are taken from independent Gaussian random variables with 
mean zero and with the variance: var (-ffoij) = 1 + % and var (Vij) = 1 -|- Sij. Note that both of 
the non-perturbed and perturbed terms have the same statistical properties. The density of states 
p{E) at the energy E is given by Wigner's semicircle law for large AT,^^) 

PiE) = 2^V4iV-£;2. (2) 

In this paper, we confine \{t) to the sinusoidal form, \{t) = Asm{ut) with A = 0.5. Eigenvalues 
of a GOE random matrix show a structure with level repulsion as a function of A. Energy spectrum 
of a system of N = 500 around the ground state, and around the center {E = 0) are depicted 
in Fig. 1 (a) and (b), respectively. Many avoided crossings are seen in Figs. 1. There are no 
degeneracies of levels although some energy levels are seen to be crossing due to the line width. 
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Fig. 1. Energy spectrum of H ioi N = 500 as a function of A around (a) the ground state and (b) the center. 



Now we consider the time evolution of a state IV') of the system; 



(3) 
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Here we take h = 1 and the initial state I'^o) is taken to be the ground state. The state IV'i) after 
a period T = 27r/a; from t = is expressed using the Floquet operator^^^ F; 

f27r/a; 



IV'i) = F|V'o), -F^Texp 



-i / n{t)dt 
Jo 



(4) 



where T means the time ordered product. Therefore the state after the nth period \iljn) is written 
as 

|V^,)=F"|V^o). (5) 

Floquet eigenvectors {\v)} form a complete set in the Hilbert space, and Floquet eigenvalues are 
aligned on a unit circle of the complex plane; 

F\i^) = e'^-\u). (6) 

We calculate F numerically by integrating the Schrodinger equation (3) for a period using the fourth 
order decomposition of time-evolution operator.^^^ The time evolution of energy is expressed in 
terms of Floquet eigenvalues and eigenstates. The energy after the nth period is expressed using 
the Floquet operator; 

ii^nmn) = (Vo|(i^^)W"|V'o) 

3. Esat as a Function of a; 

3.1 Numerical results 

We define g{i', u') and Egat as, 

g{i,,u') = {i;o\u){iy\Ho\i^'){i,'\i;o), 

E,^, = Y,g{v,v), (8) 

and rewrite eq. (7) as, 

v' ,v 

= £^sat + J2 ^0 + cc.) . (9) 

In the second term of the right-hand side of eq. (9), the expectation values of terms for {v ^ v'} 
oscillates with the period 2it/ {<i)v' — 4>v)- Thus, we see (V'nI-ffolV'n) fluctuates around £?sat- Therefore, 
Esax is regarded as the saturated value of the energy. 
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Figure 2 shows the time evolution of (^„|i7o|V'n) (zigzag curve), and the corresponding Egat 
(horizontal line) with respect to cj/tt = 0.02, 0.1, 0.2, 0.4, and 1.0 for N = 500. Figure 3 shows 
Esa,t as a function of uj obtained from five samples for N = 256, 500, and 700. We should note that 
the maximal saturated energy is zero because the energy spectrum is symmetric about zero energy 
(E=0), and this maximal energy corresponds to the high temperature limit. Here we normalize 
E'sat such that the distance between E = and the ground state is 1. Indeed in classical systems, 
the energy diffusion continues forever, which means that E'sat converges to 0. 

On the other hand, when u is very small {u = 0.027r), the time evolution of the system is 
almost adiabatic, and the average energy changes little.^^^ When uj becomes slightly larger {cj = 
O.Itt and 0.27r), transitions between levels begin to occur at avoided crossings (that is, nonadiabatic 
transitions). Thus the system absorbs energy and E{t) increases. However E{t) saturates before it 
reaches to 0. The value of £?sat gradually approaches to the center of the energy spectrum (-Egat = 0) 
as u grows larger {u = OAtt and I.Ott). These finite saturations of E'sat are caused by the quantum 
effects, and this feature is called the quantum localization.^' The quantum localization was first 
discovered in the kicked rotator model, which shows chaotic motion in the classical limit. 




pair of E{t) and E^ax denotes, from the bottom, data for u; = 0.027r, O.Itt, 0.27r, 0.4-k, and I.Ott. 

3. 2 Failure of a Picture of Independent Landau-Zener Transitions 

As shown in the previous section, the energy is saturated to a finite value. We here show 
that this saturation cannot be realized without considering the quantum interference effect. Let 
us consider the small frequency regime where transitions between energy levels take place only 
between the two levels at an avoided level crossing point. In such a region, we can introduce the 
well-known Landau-Zener picture for each transition. Wilkinson proposed a theory of the evolution 
of the energy in a random matrix model with a time-dependent perturbation.^^) His theory assumes 
that transitions take place only at avoided crossings and the transition probability is determined by 
the Landau-Zener formula. This implies that the sweeping speed of the parameter is sufficiently 
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Fig. 3. Esat as a function of uj for (a) A'^ = 256, (b) A'^ = 500, and (c) N = 700. Errorbars show the variance among 
the data obtained from five samples. 



slow, and there multiple level scatterings are not relevant. In the case that the field increases 
linearly in time, he found that this approximation of independent Landau-Zener scattering gives a 
good result. 

If we apply a simple independent Landau-Zener picture to our case , we will find, 

^sat = 0. (10) 

In what follows, we derive eq. (10). We assume that each transition occurs at an avoided level 
crossing point and the probability is given by the Landau-Zener formula. Let us consider the 
transition rate R which is the probability per unit time that the system makes a transition to one 
of the two neighboring states. This quantity R was already given for the case where the external field 
increases linearly in time.^^^ We apply the known expression to the present oscillating case. The 
transition rate R is determined by statistical distributions of the energy gaps and the asymptotic 
slopes at avoided crossing points. The average speed of the parameter X{t) is denoted by |A|. For 
the present oscillating field, |A| is roughly estimated as |A| ~ Au. In addition, let p denote the 
density of states and a denote the typical difference between the asymptotic slopes of an avoided 
crossing. Under the conditions of slow speed, M|A|/9^ ^ 1, which corresponds to the situation that 
transitions of states only occur at avoided crossings to the adjacent levels, R is expressed as, 

R CX ;iV2^1/2|3^|3/2^2 ^ ^l/2|3^|3/2^3/2_ (^^^ 

We consider the time evolution of the occupation probability for the eigenstate with energy E at 
time t, f{E,t). We can derive the diffusion equation for the occupation probability f(E,t) of the 
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state by extending Wilkinson's theory to the case that p{E) is not constant; 



5p/ ^ d (dE \ d 



dt 



dE 



where D = R/p^. Using the equation of continuity: 



dp ,T, d fdE , ^ 



eq. (12) is rewritten as, 



'dE \d\ 



df -^,dE df d 

Pm=~\^\lx'dE^dE 



""'dE 



(12) 



(13) 



(14) 



Let us consider the steady state f{E) = f{E, oo). Then the left-hand side of eq. (14) gives 
zero. Note that p is given by the semicircle law (2), p{E) = ^^V-^^ — E^ (L = 2\/iV) when N is 
very large. Therefore, f{E) is solved as 



/E i-E' 
ciexp / A;(e)de 
-L J-L 



dE' + Co, 



where cq and ci are constants of integration, and 



1 (,\,dE _^d{Dp)\ 



D 



dX 



dE J 



When we integrate eq. (15), the second term of k{e) diverges, because 



r^' 1 dDp , 
J_L Dp de J_, 



dlnp 



de — — oo, 



Dp de J_L de 

where we used D = R/p^ and eq. (11). Therefore, 

f{E) = CO. 

As a result, the average energy always becomes (the center of the energy spectrum); 

Esat = J Epf{E)dE = 0. 



(15) 



(16) 



(17) 



(18) 



(19) 



Thus we find that the idea of independent Landau-Zener scattering is not applicable at least to a 
periodic system. We cannot ignore the quantum mechanical interference effect for the quantum 
localization. 



3.3 Phenomenological Interpretation of the Quantum Localization in Analogy to the Anderson 
Localization 

We here phenomenologically interpret the quantum localization. The quantum localization 
in the present system reminds us the Anderson localization where a particle in the presence of a 
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random potential is spatially localized. The system for the Anderson localization is described 
by the Hamiltonian: 

nA = Y.Vi\i){i\+tY,\i){^\- (20) 

Here, \i) means the state of the particle at the ith site, Vi is the value of a random potential at 
the site i which is uniformly distributed with a width W, and t represents hopping to one of the 
nearest neighbor sites. In this particle system, the average hopping probability pA is expressed as 
Pa ~ The probability PA{m) that the particle exists on the mth site from the initial (0th) 

site is estimated as PA{rn) ~ ~ ^^P (— log(H^/|t|)).^^' That is, the particle is localized in the 
space exponentially. This localization is caused by the quantum interference effects. 

Let us consider the analogy to the Anderson localization in the quantum localization of the 
present model. We associate the state \i) in eq. (20) with the ith adiabatic state of Hq. Furthermore 
we focus on the frequency regime where transitions between levels take place only at avoided 
crossings. We numerically confirmed that this situation is actually realized for small oj. In this 
case, we define a transition probability pq between adjacent levels, which corresponds to pA in the 
Anderson localization, and we find the common characteristics between the present situation in the 
random matrices and the Anderson localization. That is, the localization occurs due to the quantum 
inference effect among many transitions between states. We may write the occupation probability 
(m) that a state is on the mth level as (m) ~ , where pQ is a transition probability between 
the states. Actually fast relaxation is observed numerically in pq{m) of a typical example, which 
is depicted in Fig. 4. 

0.1 , , , , , , 



qoO.05 



■l.LUiijj..i. L I.. 
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Fig. 4. The probability PQ(m) that the state is on the mth level of the system (1) for N = 500 and uj = O.Itt. It is 
obtained after 600 periods. 



If we assume an exponential form in the analogy to the Anderson localization, the saturation 
energy E'gat is written as, 

AT 

E,,t=J2PQ~'^^rn/Z, (21) 
m=l 
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where Ej^ is the mth eigenvalue of Ho and Z is the normaUzation factor, Z = 
(21) gives a finite saturation value of -Egat- 



v-AT m-l 



9 

Equation 



4. The Relevant Number of Floquet Eigenstates 

In previous sections, we showed the quantum localization in the random matrix model and 
considered an intuitive interpretation. In this section, we consider the quantum localization from 
the viewpoint of the Floquet theory. As seen in eq. (8), the quantum dynamics is determined by 
the overlaps between the initial state and Floquet eigenstates, i.e., \{i'\tpo)\'^ , and Egat is obtained 
by Floquet eigenstates. Thus we here discuss this overlap. 




01 2345678 9 101112131415 
V 

Fig. 5. The dependence of the overlap |(i^|V'o)|^ on the number v. |{j/|^o)|^ is arranged in order of the magnitude. 
The overlap decays exponentially. Here u) = O.Itt and N = 500. 

We found numerically that the distribution of Kj^IV'o)!^ decays approximately following an 
exponential function when we arrange the states in order of the magnitude of the overlap as is 
shown in Fig. 5. This property is one of the characteristics of the quantum localization. We study 
how many Floquet states are involved in the ground state, and define the minimal number iVmin of 
the Floquet states by which the initial state l-^o) is covered within a ratio r: 



1 



I i/=i 



(22) 



This quantity is the same as the quantity used in the different context by Haakc, Kus, and Scharf.^^ 
Here, we take r = 0.99. Figure 6 shows A^min as a function of lu/tt obtained from five samples for 
N = 256, 500, and 700. 

We introduce the fitting function for N^ini^^) in the form: 



iVmin(a;) = iV*i, 1 



1 



(23) 



[(a;/u;*)« + l]''J 

This function fits the data quite well. In eq. (23), N^^^, a, b, and oo* are determined from the least- 
squares method. The results are shown in Table I. From eq. (23), N^i-a{Lo) is extremely small when 
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Fig. 6. A'^min as a function of w/tt for (a) N = 256, (b) A'' = 500, and (c) A'' = 700. Errorbars show the variance 
among the data obtained from five samples. The fitted fine and also linear fitting in the intermediate region are 
shown in each graph. 

a; <C o;*. This is partially because the energy gaps in lower energy region are so large comparing 
with uj that transition to higher levels from the ground state seldom occurs (see Fig. 1(a)), and 
therefore the system behaves almost adiabatically.^^^ On the other hand, N-amiit^) is well fitted by 
a linear function (dashed line) for the region of a; ~ a;* . 

So far, we observed localization in the energy space and also in the space of Floquet states. 
Now we discuss on localization length of the two localizations. In particular, we focus on the linear 
dependence of A?inin(t^) on uj in the small u region. We employ a phenomenological argument 
using the analogy to the Anderson localization as explained in §3.3. Let us express the transition 
probability pQ using some function r(a;) in the form: 

PQ~exp[-r-i(u;)] . (24) 

The occupation probability PQ{m) that a stationary state is on the mth level is expressed as 
pq{rn) ^ Pq ^ exp [— r~^(a;)m] . Thus we can regard T{oj) as the localization length in the energy 
space. Now let / denote the localization length in the Floquet space, i.e., |(z/|'!/'o)P — ce~'^/', where 
c is the normalization constant. From the definition (22), iVmin is expressed using the length I as 

/ 



A^min --J^ -r)<X.l. 

Since we now focus on the linear dependence of A^min on w, we write, 

-Amin = au + (3, 

where a and /3 are constants. We note that /3 is negligibly small. Then we have 

I oc u. 



(25) 



(26) 



(27) 
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Table I. Values of the parameters in eq. (23) obtained by the least-squares method. 
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variables 


dim(256) 


dim(500) 


dim(700) 




0.79 


0.80 


0.80 


a 


2.3 


2.6 


3.0 


b 


0.67 


0.47 


0.54 




0.15 


0.10 


0.11 



As has been assumed in the the kicked rotator model to discuss the diffusion constant,^^^ we assume 
the proportionality between I and r(u;), 

I oc r(a;). (28) 

Then, we obtain from eqs. (27) and (28), 

r{uj)=co/h. (29) 

Here h is an unknown amplitude which may depend on N, A, the typical size e of gaps at avoided 
crossings, and the difference a of the two asymptotic slopes at an avoided crossing. This u depen- 
dence is consistent with that of the Landau-Zener formula, because eqs. (24) and (29) give 

(30) 

Actually, we numerically confirmed that transitions between levels take place only at avoided 
crossings in the region where A'min linearly depends on oo. Thus we suppose that the amplitude 
of local transition probability is originated in the Landau-Zener transition, although the phase 
interference has an important effect on the global phenomena. 

5. Summary and Discussion 

We have studied numerically the time evolution of the average energy of the GOE random 
matrix model under a perturbation of a sinusoidal function of time. In §3, we showed that the 
quantum localization occurs in this model. We have also studied the frequency dependence of the 



PQ ~ exp 



h 
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saturated energy Esat, and we showed that we can not rely on the picture of independent occurrence 
of Landau-Zener transitions in the present case because the suppression of the quantum diffusion 
is essentially due to the phase interference effect of quantum mechanics. We discuss this quantum 
localization in analogy to the Anderson localization. In §4, we introduced the relevant number of 
Floquet states iVmin, and proposed a form of Ni^m as a function of u. Ntain is almost zero when 
u is nearly zero, and it saturates when to is sufficiently large. In the intermediate region, Nja.m 
grows linearly as u increases. This linear dependence of iVmin on u implies that the localization 
length I in the Floquet eigenstates also depends linearly on u (eq. (25)). On condition that F oc Z, 
the localization length F of the eigenstates of the Hamiltonian is also proportional to u. This 
dependence of F on a; implies the Landau-Zener mechanism governs the local transitions. 
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